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Abstract 

Using the spherically symmetric self-consistent Green’s function method, we consider 
thermodynamic properties of the S = 1/2 J\-J^ Heisenberg model on the 2D square lattice. 
We calculate the temperature dependence of the spin-spin correlation functions c r = (SqS*), 
the gaps in the spin excitation spectrum, the energy E and the heat capacity Cy for the 
whole J\ .^-circle, i.e. for arbitrary 92 , J\ = cos(tp), J 2 = sin(<p). Due to low dimension there 
is no long-range order at T 7 ^ 0, but the short-range holds the memory of the parent zero- 
temperature ordered phase (antiferromagnetic, stripe or ferromagnetic). E((p) and Cy(<p) 
demonstrate extrema ’’above” the long-range ordered phases and in the regions of rapid 
short-range rearranging. Tracts of c r ((p) lines have several nodes leading to nonmonotonic 
c r (T) dependence. For any fixed <p the heat capacity Cy(T) always has maximum, tending 
to zero at T — >• 0, in the narrow vicinity of ip = 155° it exhibits an additional frustration- 
induced low-temperature maximum. We have also found the nonmonotonic behaviour of the 
spin gaps at (p = 270° ± 0 and exponentially small antiferromagnetic gap up to (T < 0.5) for 
<P > 270°. 

PACS codes 

75.10. Jm Quantized spin models, including quantum spin frustration 

75.10. Kt Quantum spin liquids, valence bond phases and related phenomena 

75.30.Kz Magnetic phase boundaries 
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1 Introduction 

Investigation of low-dimensional quantum magnets have attracted much attention during the 
last years (see [1] for a review). Frustrated two-dimensional (2D) and quasi-2D magnets are of 
particular interest, as they demonstrate strong quantum fluctuations effects. The 2D spin-1/2 
J 4 -J 2 quantum Heisenberg model is a conventional tool for the investigation of frustration effects 
and quantum phase transitions (see, e.g., m)- 

Cuprates and numerous other quasi-2D compound with antiferromagnetic (AFM) nearest- 
neighbour (NN) and next-nearest neighbour (NNN) couplings J\ > 0, J 2 > 0 have been investi¬ 
gated experimentally for years mm 

This class of systems has been recently expanded by several magnetic materials with a fer¬ 
romagnetic (FM) NN coupling J\ < 0 and a frustrating AFM NNN coupling J 2 > 0, e.g., 
Pb 2 V0(P0 4 ) 2 [IMS], (CuCl)LaNb 2 0 7 HE], SrZnV0(P0 4 ) 2 USES EE], and BaCdV0(P0 4 ) 2 
ns im mg. The frustrating J 2 is believed to be large enough to drive these materials out of 
the FM phase. There are also materials, such as K 2 CuF 4 , Cs 2 CuF 4 , Cs 2 AgF 4 , La 2 BaCu05, and 
Rb 2 CrCl 4 , pm [TT1 with insufficiently strong frustrating AFM NNN interaction. 

The general picture can be seen from Figure El where the phase diagram of 2D Ji-J 2 model 
in the classical limit S —> oo is complemented by the positions of several experimental systems. 
Hereinafter the Ji-J 2 -circle is defined by J\ = cos(ip), J 2 = sin(ip). In the classical limit only 
three phases are realized — all with long-range order (LRO) — AFM, FM and stripe (in the 
quantum case S' = 1/2 disordered phases appear between stripe and FM, as well as between stripe 
and AFM). 

The Hamiltonian of the model has the form 

H = J, Y, SiSj + h Y, SiSj (1) 

<D> [D] 



where (Si ) 2 = 3/4, (i, j) denotes NN bonds and [i, j] denotes NNN bonds of the square lattice sites 

i,j- 

The theoretical investigation of the model in the first quadrant J 4 > 0, J 2 > 0 was detonated by 
the HTSC breakthrough and led to innumerable number of papers (see, e.g. [241125] and references 
therein). In a nutshell, the generally accepted result is the following. At T = 0 the system 
undergoes two successive phase transitions: from AFM LRO to disordered phase and then to 
stripe LRO (see, for example, recent calculations [26] and references therein). The nature of these 
quantum phase transitions and the detailed structure of the disordered state remains debatable. 

The unfrustrated FM case (J 4 < 0, J 2 = 0) has been also widely investigated, e.g., by the 
modified spin-wave theory [271128] . renormalization group approaches [291130] , the quantum Monte 
Carlo method [3T1134] and by a spherically symmetric self-consistent approach (SSSA) — spin- 
rotation-invariant Green’s function method (RGM) in alternative notation — [35H22]- 
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Figure 1: (Color online) Circle phase diagram of the 2D J1-J2 Heisenberg model in classical 
limit; points on the circle correspond to exchange parameters of several layered compounds, data 
from [HH] 


Recent experiments with vanadates and related compounds have stimulated theoretical studies 
of the Ji-J 2 model with < 0 and frustrating J 2 > 0. [4QH5Tj . ft was found that in the second 
quadrant there is also a disordered ground state between FM and stripe. Rough estimates for the 
transition points in both quadrants are J 2 ~ ±0.4Ji (AFM —y disordered, FM —>■ disordered) and 
J 2 ~ ±0.7 A (disordered —> stripe ). Note that for the classical model (S —$■ 00 ) there are only 
two transitions at points J 2 = ±0.5 A (AFM —> stripe, FM —$■ stripe). 

So, one has several experimental points settled on the upper half of the A _ A-circle and a set 
of theoretical methods, each being tuned for a particular parameter region. A unified approach, 
which can describe the entire picture, both for the ground state and the thermodynamics of the 
model, is obviously desirable. It is also interesting to look at the lower half of the circle (A < 0), 
though still experimentally unobtainable. 

The SSSA proved to be the appropriate approach. SSSA preserves the spin SU(2) and trans¬ 
lation symmetries of the Hamiltonian and allows: 

i. to satisfy automatically the Marshall and MerminWagner theorems 

ii. to describe at T = 0 (when the LRO is possible) the system states both with and without 
LRO in the frames of one and the same approach 

iii. to find the microscopic characteristics such as the spin-excitation spectrum cu(q), the spin- 
gaps T-dependence and the explicit form of the dynamic susceptibility y(q, u>, T); to go beyond the 
mean-field approximation by introducing damping in the expression for the spin Greens function 

G(q,cu) [Zj. 

Note, that for the ID spin-1/2 Heisenberg ferromagnet it was shown that the RGM reproduces 
Bethe-ansatz results [371152] , 

Let us also mention that SSSA has been applied to another lattice geometry, S' > 1/2 [53], the 
systems with the anisotropic spin exchange — iron pnictides AFe 2 As 2 (A = Ca, Sr, Ba) — [ 54] 
and to 2D A _ A _ A model [711551136] . 
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We have recently considered the ground state of the model dTJ) Ji—J 2 circle [57]. In particular, it 
was shown that at T = 0 the transitions between all ordered and disordered phases are continuous, 
except the transition FM—^AFM at J\ = 0, J 2 = —1. 

In the present work we turn to the nonzero temperatures and consider thermodynamic prop¬ 
erties. We calculate the temperature dependence of the spin-spin correlation functions, the gaps 
in the spin excitation spectrum and the heat capacity. As in [57], the model is treated in SSSA. 

The paper is organized as follows: in Sec. [2] and Sec. [3] we briefly remind the spherically 
symmetric self-consistent approach for two-time retarded Green’s functions and the results of ra 
for T — 0. In Sec. [4] we present the results for the thermodynamic properties. Discussion and 
summary are given in Sec. [5] 

2 Spherically symmetric self-consistent approach 

As already mentioned, the calculations are performed in SSSA for the spin-spin Green’s function 

EESlEBlSSlSSlEra. 


G„ m = (S'lSi.Uii = ~ijdte**({Si(t), SJJ) (2) 

0 

Due to the spherical symmetry, only the Green’s functions diagonal with respect to a = x, y, z 
are nonzero, mean cite spin is zero ( G zz = G xx = G yv \ (S@) = 0, /? = x,y,z). There are 
three branches of spin excitations degenerate with respect to /3. Because ([S^, S^J) = 0 and 
[Sn, Sm] = ^nm£a/3 7 S'n, We llaVe for G nm 


wG nm = {[Si, SJJ) + ([Si, ff]|Si,}„ +u = i ^ J b ((Si +b Si - Si +b Si)|Si,}„ +u ( 3 ) 

b=g,d 

where J g — J \, Jd = T 2 . The second differentiation step leads to three-site Green’s functions: 


™ V J b {(Si +b Si - Si +b Si)|Si,). +i , = 2 V J b C b (i„ +b , m - 5 „. m )- 


E J b(G„ + b. m - GJ - V ./ b ./ b .[(Sf +b Si_ b ,S„ s |Si,)„ + „- 


b=g,d 

20 Q0 


b,b'#z;b+byO 
20 


- (SZ +b Sl h 'S z n \S z m ) u+u + (^ +b ^ +b+b ,^|^) w+i , - (^b^b+b'^L+J} (4) 

where c b = {S* +b S$. 

In the mean-field approximation, the subsequent procedure amounts to decoupling the chain of 
the equations of motion at the second step using the triple-site term splitting of the characteristic 
form 

‘- > n+g 1 +g 2 ‘^n+g, sz « ( , 5'n+g 1 +g2^'n+g 1 

+^7 (S l n+Si S^\ *S'n+ gl +g 2 ) + a gl+g2^J7 (^n+ Sl + S 2 ^n) ^n+s 1 


(5) 


where a g and a gl + g2 are vertex corrections. In the general case vertex corrections can depend 
on cite indices, hereafter we use the simplest one-vertex approximation [HH11HB10BI1SZ], he. all the 
vertices are taken to be equal. We emphasize that in the case of S = 1/2 the splitting procedure 
for each term is unambiguous because in (|41) /j A z, and the average involves two spins with (3 ^ z. 
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After the Fourier transformation 


^ = 7^ £ e-i,r ^ M 

solution of equations (j3J) and (J5J) leads to the Green’s function G(q, oj) = (5q|<S£ q )u> = — x(q, cc) 

G(q,«) = , /" ,, (7) 


In the case of J\ 


(^ 2 - ^ 2 ) 

-J 2 model the expressions for A q and cUq are the following: 

F q = -8 [JiC 5 (l - 7 g) + J 2 c d (l - 7 d )]; 

= 2[(7iAd + 7 2 A" 2 ) - (7 3 AT 3 + 74A4) - (75JF5 + 7 6 AT 6 )] 

The variables K\... A£ involved in the expression for the spectrum are given by 

K\ = J\J 2 K g d + 12 J\c g + 1 + K gg ] K 2 = J\J 2 K g d + 12 . 7^7 + 1 + Kdd 


16 J 2 c s ; 

K 4 = 16J£c d ] 

— 16JiJ 2 c ff ; 

IQ = 

k 99 = 

^ ) C r , Kdd 

^ ) C r , Kgd ~ 

- £ 


r =gi+g2; 

r=di+d 2 ; 

r=g+d 


gl^”g2 

di^—d2 



( 8 ) 

(9) 

( 10 ) 

( 11 ) 

( 12 ) 


71 = 1 -7»; 72 = l-7d; 73 = 1 -7gi 74 = 1 -7di 75 = (l-7s)7d5 (13) 

7 e = (1 - 7 d ) 7 5 ; 7 5 = -(cos q x + cos q y ); 7 d = cos q x cos <37 (14) 

where c r stands for the correlators, renormalized by the vertex correction c r = ac r (in (J9]) we 
rearranged the contributions to the spectrum in a different manner, than in [ IF . TTTl). 

The Green’s function G(q, oj) involves the correlators c r for the first five coordination spheres 
r = g, d, 2g, g + d, 2d, which must be evaluated self-consistently in terms of G(q, a;). In addition, 
G(q, co) must satisfy the spin constraint c r= o = {S?S?) = 1/4 (sum rule). These conditions are 

c r = — c q e iqr ; Cq = (R q £* q ) = — / duj ( 2m ( ) + 1) lrnG zz (cc, q) ; (15) 

n 0 


c r=0 = 1/4 


1 1 
7T iV 



du (2m(w) + 1) lmG zz (cc, q); 


(16) 


The system of self-consistent equations is then solved numerically. Hereafter all the energy- 
related parameters are set in the units of J = \J Jf + Jf. 

In the general case, in the framework of (I3j), (HI), both short-range order (SRO) and long-range 
order (LRO) states can be realized. Because the dimension is equal to two, only SRO are possible 
at T 7 ^ 0, and both possibilities can take place as T —> 0 (LRO is characterized by nonzero 
spin-spin correlators at infinity (S^Sq ) r ^. 00 ). 
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Figure 2: (Color online) Condensate c con d (black bold line) and correlation functions on the first 
three coordination spheres as functions of the angle parameter ip = cos ip, J 2 = sin ip) — c g 
(blue), Cd (green) and c 2g (red); bold lines correspond to the zero temperature, solid lines — to 
T = 0.3, 0.4, 0.5, dashed lines — to T — 0.6, 0.7 and dotted lines — to T — 0.8, 0.9. 

3 The ground state properties 

For T 7 ^ 0 the system is always in the SRO state (the spin-liquid state). For small T, the function 
c q is always peaked near the point qo where the spin gap is minimum (except for the zero point 
q = 0 , where spin gap is always zero, but c q is not peaked because F q=0 = 0 ). 

Two cases are possible as T —y 0. In the first case, the system remains in the SRO state, the 
gap is not closed, and the contribution of m(co) to c q in (fl5l) vanishes as T —» 0. 

In the second case, as the temperature decreases, the system passes into the LRO state, and 
the gap goes to zero at a point q 0 , which for J\ — J 2 model can be either qo = Q = ( 7 r, 7 r), 
(corresponding to antiferomagnetic LRO) or q 0 = X = ( 7 r, 0 ); (0, 7 t) (stripe LRO). The zero 
gap leads to the existence of condensate c con d, which determines spin-spin correlation function at 
infinity. That is nonzero condensate means the LRO existence at zero temperature. 

The case of FM LRO is somewhat different. The condensate at the point q 0 = T = (0, 0) 
does appear for T —* 0, but the spin gap at q 0 is closed for any temperature. In this case the 
transition to LRO is governed by the spectrum transformation near qo — from linear to quadratic 
(see detailes in ra>- 

Note that the third exchange interaction J 3 being added to the model changes the LRO picture 
qualitatively — the helical LRO becomes possible. In the J\ — J 2 — J 3 model the condensate 
peak point in the structure factor can be located not only at T, Q, or X, but also at arbitrary 
incommensurate point on the side or diagonal of the Brillouin zone [7i [46ll55lI 56 [ [60] . 

The resulting picture for the whole region of J\ and J 2 exchanges (” J\-J 2 circle”) EH is shown 
in Figures I2ll4l As seen in the figures, four types of the ground state are possible: AFM, stripe, 
FM and the disordered spin-liquid state (two different areas — SL 1 and SL 2 ). Note, that in the 
frames of SSSA the transition SL 2 —>RM appears to be continuous (W)- 
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Let us remind, that in SSSA mean spin projection is always zero ( S £) = 0 and any possible 
ground state preserves the whole — translational and spin SU(2) — symmetry of the Hamiltonian. 
The LRO, if it does exist (this possibility is open — LRO is possible only at T = 0) is determined 
by the spin-spin correlation function at infinity. 

At ip = 0 (Ji = 1, J 2 = 0), as it the classical limit, the ground state is AFM. Quantum 
fluctuations destroy LRO with the increasing ip (that is increasing J 2 ) and the system transforms 
to the disordered state SL 1 . Note, that different alternative states are considered to be competitive 
in the locality of LRO disappearance at T = 0. These are in particular columnar and box phases 
which preserve the SU( 2 ) symmetry, but brake the translational one (see [ 6 T] for recent review). 
The disordered state in SSSA is always spin-liquid-like in the above noted sense and it can not be 
distinguished from the mentioned alternatives with nearby energies. 




c 

Figure 3: (Color online) Polar diagrams for the absolute values of the first (a), second (b) and third 
(c) correlation functions (|c 9 |, |q|, and |c 2ff | correspondingly); outer lines correspond to T = 0, 
inner lines — to T = 0.3 -t- 0.9. 


In the vicinity of p = 7 t/2 there exists another ordered state, with stripe LRO type. Then 
with the increasing <p the spin liquid appears again (SL 2 , its SRO differing from that of SL 1 ). 
After SL 2 the system continuously transforms to FM-LRO state. And the last transition — at 
ip = 37 t /4 restores the AFM LRO. In the frames of SSSA all the mentioned transitions at T = 0 
are continuous, except the last one. 

4 Thermodynamic properties 

The results of our calculations for the spin-spin correlation hnctions on the first, second and third 
coordination spheres (c g , Cd and c 2g correspondingly) at different temperatures (0 < T < 0.8) are 
presented in Figure [2] 
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Figure 4: (Color online) Spin gaps Aq (blue) and Ax (green) at points Q = ( 7 r, 7 r) and X = 
(7r,0), (0,7r) as functions of the angle parameter tp (Ji = cos <p, J 2 = sin ip)] bold lines correspond 
to the zero temperature, solid lines — to T — 0.3, 0.4, 0.5, dashed lines — to T — 0.6, 0.7 and 
dotted lines — to T = 0.8, 0.9. Black bold line corresponds to the value of condensate function 
Ccond (hi altered scale). 


As it is seen from the figure, at nonzero temperatures the SRO of the disordered state holds 
the memory of the parent zero-temperature ordered phase. The region overlying the AFM phase 
is characterized by AFM-like correlators c g < 0, \c g \ > Cd> C 2 g > 0, the region above stripe phase 
- by stripe-like ones c ( i < 0, C 2 g > 0, \cd\ > C 2 g > \c g \, and the area above FM phase - by FM-like 
c g > Cd > C 2 g > 0. In the intermediate regions SRO transforms from one limit to another. 

It is natural to expect simultaneous decrease of all the correlators absolute values with the 
temperature growth. In the main so it is. There are, however, two important exceptions with 
inverse temperature dependence. It is firstly the nearest correlator c g in the region of small J\ - 
above the stripe phase, in part of SL 1 and in the whole SL 2 region; secondly, NN correlator Cd in 
narrow vicinities of the nodes - in the middle of SL 1 and near the transition SL 2 —>FM (see Sec. [5] 
for the explanations). 

The general view of the correlators temperature behaviour can be also seen in Figure Eh, 
Figure Eh and Figure Eh, where polar diagrams for the correlators absolute values are shown. 
These figures allow to compare the SRO structure with the initial classical phase diagram. In 
particular, it is clearly seen, that in the region corresponding to classical stripe phase the nearest 
correlator \c g \ is an order of magnitude smaller than the next-nearest one |q|• The regions of c g 
and Cd inverse temperature dependence can be also seen. 

The gap in the spin excitations spectrum at zero point T =(0,0) of the Brillouin zone is zero 
Ar = 0 in any phase at any temperature. The temperature dependence of spin gaps Aq and 
Ax at two another symmetrical points — AFM point Q = ( 7 r, 7r) and (equivalent) stripe points 
X = (0,7r), (7T, 0) — is shown in Fig. [7] (the corresponding polar diagrams - Fig. Eh Fig. Eh). 

Figure El shows the results for heat capacity Cy{p>) at different temperatures. Two next figures 
- Figures [7] El — show the detailed temperature dependencies for several values of ip (in the 
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Figure 5: (Color online) Polar diagrams for the spin gaps Aq (a) and Ax (b) (compare with 

FigHD- 


temperature intervals, where the adopted computational scheme leads to satisfactory convergence). 
The two mentioned figures demonstrate in particular, that the temperature curve for heat capacity 
has a peak for any chosen value of cp. 

5 Discussion and summary 

Hereinafter we discuss the most interesting properties of the phase diagram obtained. Let us 
remind, that the LRO is absent at T ^ 0 and all the ” circle” 0 < cp < 2tt is covered by spin-liquid. 
Nevertheless we will classify the 0 phase diagram regions by the LRO type at T = 0. 

5.1 AFM order 

5.1.1 

A dramatic difference of the T — 0 LRO strength (the value of spin-spin correlator at infinity) at 
points cp — 0 (Ji = 1, J 2 = 0) and (p = 37t/2 + 0 (J\ — +0, J 2 = — 1) is worth noting. 

That is the infinitesimal AFM NN exchange with the FM NNN (diagonal) exchange equal 
to unity leads to much more strong LRO than in the case of conventional AFM with NN AFM 
exchange equal to unity. The SRO difference between these two points (both at zero and nonzero 
temperatures) is not so large. 

5.1.2 

In the interval between cp = 3 tt/2 + 0 and <p ~ 330° (Figure 0] and Figured) the AFM spin gap 
is small (Aq -C T) up to high enough temperatures (T < 0.5). At the conventional AFM point 
cp = 0 the gap becomes exponentially small at much lower temperatures (T < 0.05 <C 0.5). When 
Aq <C T the gap Aq determines the AFM correlation length £afm ~ Aq 1 . So the correlation 
length increases greatly from the point cp — 0 to the interval 37 t/2 < cp < 330° (that is consistent 
with the mentioned LRO evolution). 

Note, that, as it can be shown, at cp = 0 the exponentially small gap at high temperatures 
T ~ 0.5 is realized for the spin S > 1. 
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5.2 Spin liquid 

5.2.1 

Let us note the dissimilarity of the spin liquid evolution (with growing <p ) in two areas, corre¬ 
sponding to zero-temperature regions SL 1 and SL 2 . The SL 1 liquid at T = 0 exists in the angle 
range 0.051 < ip < 1.111 (2.92° < ip < 57.30°), it is between AFM and stripe phases. The SRO 
in SL 1 smoothly transforms from the left bank to the right one, from the AFM-type (c g < 0, 
\c g \ > Cd > C 2 g > 0) to the stripe-type (cd < 0, C 2 g > 0, \cd\ > C 2 g > \c g \). That is not the case 
for SL 2 liquid (2.141 < ip < 2.712 (122.67° < ip < 155.39°) at T = 0, the region between stripe 
and FM phases). Almost in the whole SL 2 area the SRO is stripe-like, in particular, Cd remains 
negative. The absolute value of q almost everywhere, except tiny region near ip = 2.712, is larger 
than the nearest neighbour correlator c g . And only close to FM border the drastic restructuring 
of the correlators takes place. At T ^ 0, as it was repeatedly mentioned, there is no LRO, but all 
the above statements concerning SRO do hold. 

5.2.2 

As it was noted earlier, there are regions of the phase diagram with anomalous temperature 
behaviour of the correlators at fixed ip (absolute value growing with temperature or nonmonotonic 
behaviour). 

This anomaly for c g correlator realizes in the region from the middle of SL 1 phase through the 
stripe-phase and all the SL 2 up to the transition to FM. The anomalous regions for q correlator 
- narrow areas near the nodes — are in the middle of SL 1 and near SL 2 —>RM transition. The 
reason for Cd anomaly is obviously the following. For different temperatures Cd(<p) changes sign at 
different points (p. The normal temperature behaviour of the correlator would pass to the normal 
one after crossing <p axis only if the node of the correlators cone should be exactly on the ip axis 
(that is Cd(<Po,T) = 0 for any T). But generally there is no physical reason for this statement. 

The reasons for other mentioned temperature anomaly — the c g (<p,T) behaviour — are not 
so obvious. Presumably it is connected to the rapid SRO rearrangement. 

5.3 Stripe order 

The most interesting point of the area is <p = 7t/2 {J\ = 0, J 2 = 1). At this point the lattice is 
decoupled into two non-interacting sublattices. One can see from Figure [21 that, as it should be, 
at any temperature Cd(j r/2) = c 9 (0), C 2 g (tt/ 2) = Cd(0). The decoupling means that c g {ir/2,T) 
is strictly zero. That is why the cone of c g correlators with anomalous T-behaviour retains the 
anomalous behaviour after crossing the node on ip axis. 

At the same point the AFM gap Aq(y? = 7 t/ 2) = 0 for any temperature (see Figure [I]), though 
AFM LRO at any temperatures, including T = 0, is absent. 

Formally it follows from the analytical expression for the AFM gap Aq (ip = vr/2) ~ J\. 
The naive explanation is that the system is degenerate with respect to mutual rotation of the 
sublattices, that is the transfer of spin excitation to the neighbouring cite costs no energy. 

All the same is true for the point ip = 3 tt/2 (J\ = 0, J 2 = —1). Though is this case there is 
AFM LRO at p> = 3n/2 + 0 and T = 0 (and FM LRO at ip = 3it/2 — 0, T = 0). 
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5.4 FM order 

5.4.1 SL 2 —>■ FM transition 


As it was noted earlier, it was shown in m, that SL 2 —>■ FM transition at T = 0 is continuous, 
though it takes place in the very narrow p interval. At T ^ 0, as it is seen in Figures dd hi the 
vicinity of this transition the c^g temperature dependence is nonmonotonic and the temperature 
dependencies of other correlators and the gaps are inverse. The heat capacity (as function of ip) 
at any temperature has sharp minimum near this transition (Figures [6HH]). 

5.4.2 FM —> AFM transition 

This transition takes place at p = 37t/2. At this point the lattice is splitted into two noninteracting 
sublattices. At (p — > 3n/2 — 0 there is no frustration with respect to the FM order, at (p —> 37t/2 + 0 
- no frustration with respect to the AFM order. Therefore it is physically obvious that in the 
quantum limit at T = 0 a transition between these two phases is of the first order and there is no 
spin-liquid area between FM an AFM. Our previous calculations m do confirm this evidence. At 
nonzero temperatures the transition is obviously continuous. But as it is seen from Figures dd 
the correlator c g (p ) rapidly transforms and changes the sign. Other correlators and the spin gaps 
demonstrate nonmonotonic dependence on ip , the heat capacity has sharp minimum. 

Note once more, that near ip = 270 + A p the AFM gap Aq is exponentially small up to 
T ~ 0.5. Aq 10 -4 - 1(T 5 for Ap ~ 5°. That is why our calculation can not reproduce the fine 
structure of the correlators and the heat capacity for ip = 270 + A p at low temperatures. 

5.5 Specific heat 

As it is seen from Figures dd the heat capacity for any p tends to zero at T —» 0. The reason is 
the stabilization of all the correlators (and the energy) at low temperatures. At any fixed p the 
heat capacity T-dependence, as it should, has the maximum, varying in its height and position 
(see Figure [Tf Figure IHJ). 

It is also seen from Figured that there always exists a local minimum of the function Cy{p) 
inside the areas above the ordered phases (AFM, stripe and FM), if T is not too high. The reason 
is that the correlators (and the energy) in the corresponding areas weakly depend on p. These 
minima correspond to local minima of the energy E(p). 

Three other local minima of Cy(p) are observed, on the contrary, in the regions of neighbouring 
short-range orders rivalry, where the correlators are rapidly rearranging, and the energy also weakly 
depends on p. These minima correspond to local maxima of the energy E{p). 

Five of the mentioned points of Cy(p) local minima are marked by diamonds in Figure d The 
heat capacity temperature dependencies C V (T) at these points and symmetrically neighbouring 
points are shown in Figures d d 

Figure [7] corresponds to p in the middle of ordered (at T = 0) Stripe and AFM phases. It is 
interesting that the heat capacity Cy(p,T) is symmetric on p relative to the minima of Cy(p) 
and E(p) in a wide T range. The situation in the region above the FM phase is slightly different, 
because, as it is seen from Figured Cy(p) moves with temperature. 

Figure d corresponds to p in the regions of the strong frustration where different short-range 
orders concur. The heat capacity Cy(p,T) is symmetric on p relative to the minimum of Cy(p) 
and maximum E(p) in a wide T range. The case of p = 155° ± 5° is obviously notable. As in 
other cases symmetrical satellite Cy(T) lines are similar, but the central line Cy(p = 155°, T) 
demonstrates additional maximum at low T . This maximum was found and discussed in [48]. The 
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Figure 6: (Color online) Specific heat Cy (blue) and energy per site E (green) as functions of 
the angle parameter ip = cos ip, J 2 = sin <p). Bold green line corresponds to the energy of the 
ground state (T = 0). Solid lines correspond to T = 0.3, 0.4, 0.5, dashed lines — to T = 0.6, 0.7 
and dotted lines — to T = 0.8, 0.9. Black bold line corresponds to the value of condensate function 
Ccond (in altered scale). The Cy(T ) slices for p> values marked by diamonds and tics are shown in 
Figures EHHl 



Figure 7: (Color online) Families of the specific heat Cy(T) curves near values of angle parameter 
ip (Ji = cos 99 , J 2 = siny?) corresponding to Cy((p) minima above FM and AFM phases (see 
Figure [6]). 
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Figure 8 : (Color online) Families of the specific heat Cy{T) curves near values of angle parameter 
ip (Ji = cosip, J 2 = sin ip) corresponding to Cy(ip) minima above the regions of SRO rearranging 
(see Figure [ 6 ]). Cy(T,ip = 155°) demonstrates two maxima - see text. 


low-temperature frustration-induced maximum Cy{T) was also found in ID case for S = 1/2 and 

s =1 m- 


5.6 Fine tuning of the method 


As it was noted above, all the calculations in the present work has been carried out in the straight 
and simple approximation, without any tuning parameters. The tuning in SSSA is commonly made 
via different manipulations with vertex corrections [3EIlEaESlSHlSniE2lElMlE2] or accounting 
for complex structure of the Green’s function, that is, considering the polarization operator (in 
particular, accounting for damping of spin excitations) [71 [551156] . 

All such complications obviously affect the results. At very low T the obtained quantitative 
differences can amount significant values. Figure E] demonstrates, that even the simplest self- 
consistent accounting for the damping in the frames of the Green’s function 


G zz 


(u, q) = 


UT — Wq + iu'y 


(17) 


(7 — damping parameter) shifts the borders of the disordered phase between AFM and stipe 
phases at T = 0. 

Nevertheless, our estimates and comparison of the available data confirm, that the reasonable 
degree of tuning does not lead to any topological modifications of the phase diagram. 


5.7 Summary 

To summarize, in the present work thermodynamic properties of of the 2 D J\ — J 2 S = 1/2 Heisen¬ 
berg model are considered for the entire phase diagram the frames of one and the same approach 
- spherically symmetric self-consistent approach for two-time retarded Green’s functions. 
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Figure 9: (Color online) Effect of the damping parameter 7 on the phase boundaries of the spin 
liquid SL 1 . 

This work is supported by Russian Foundation for Basic Research, grant 13-02-00909a. 
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